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Abstract 

CN) Properties of systems with majority voting rules have been exhaustingly studied. In this work 

we focus on the randomized case - where the system is initialized by randomized initial set of seeds. 

i—5 Our main aim is to give an asymptotic estimate for sampling probability, such that the initial set 

CS1 of seeds is (is not) a dynamic monopoly almost surely. After presenting some trivial examples, we 
present exhaustive results for toroidal mesh and random 4-regular graph under simple majority 

i — i scenario. 



1 Introduction 



Idea of majority voting is commonly used to resolve many problems related to achieving consis- 
tency between different parts of distributed computation. For example, majority voting is used to 
preserve data consistency when updating copies of the same data. Also, it is quite common to use 
majority voting to resolve inconsistencies in distributed database management. Majority based 
systems were also successfully used by Agur to examine the plasticity and precision of the immune 
response in [17]. 

The model for the system is as follows: Let G = (V, E) be simple undirected graph of size 
n. Every vertex has its color which is either black or white, and represents the state of the 
node (for example black = faulty, white = not faulty). By S we shall denote the set of black 
vertices at the beginning of the process. These vertices are also called seeds. By evolution of 
such system (or, the coloring process) we shall mean the synchronous process where in each step 
each vertex adjusts its color according to colors of its neighbors and its internal contamination 
and decontamination rules. (De) contamination rules determine under what configuration of its 
neighbors the white (black) vertex turns black (white). In this work we focus on the case when 
there is no decontamination rule, and the contamination rule is the simple majority rule. This 
means that the white vertex turns black if at least half of its neighbors are black and there is no 
possibility for black vertex to turn white. Set of seeds S is called dynamic monopoly, (or dynamo 
for short) if the corresponding coloring process leads to monochromatic black graph. For more 
rigorous definition, see section 2.1. 

A significant attention was given to this model, and many interesting results were obtained. 
Probably the most basic (but certainly not trivial) question asks for determining the minimal 
cardinality of a dynamo on a fixed graph G [4, 5]. Another interesting parameter of a dynamo 
besides its cardinality is the time that is needed for contamination to spread. This was analyzed 
for the first time in the literature in [4, 5]. Several works are related to more advanced topics such 



as decontamination of the system by external agents [7]. Finally, [9] defines the terminology of 
immune subgraphs and asks how the immune subgraph of a certain graph looks like. 

All these tasks were solved for small class of graphs. Close attention is given to the ring and 
its several modifications [4], as well as to the torus and its modifications [5] (toroidal mesh, torus 
cordalis, torus serpentinus). Many results exist also for hypercube and binary tree. 

Other possible questions ask about the minimal cardinality of a dynamo on arbitrary graphs. 
It is proved that on general directed graph the minimal dynamo has at most . 72 7 1 V^ | vertices. It is 
likely that this result will be improved. On the other hand, in the case of general undirected graphs 
the minimal dynamo consists of at most |_V/2J (+1) vertices (depending on the exact specification 
of contaminating rules), this result is proved to be sharp [18]. 

There exist some interesting results about such systems even in the "most general" case, where 
contamination rule is that the vertex is contaminated if at least fraction a of its neighbors are 
black [20, 21, 22]. These works ask under what condition at least some fraction 5 of all vertices is 
turned black. 

Finally, we mention several works that we find most related to the paper. For a randomly 
chosen set of initial black vertices, Gleeson and Cahalane [21] gave an exact formula for the 
expected fraction of black vertices at the end in tree-like graphs. In [20, 23] authors gave their 
estimate for minimal number of black vertices needed for re-coloring of at least fraction 8 of all 
vertices on Erdos-Renyi random graph. For more thorough survey, see [10, 8]. 

Motivation for studying coloring process induced by random set of seeds is quite straight- 
forward. For example, considering vertices as computing nodes, each node can fail with probability 
p, independent on failure of other nodes. Although this looks like a very common scenario, there 
exists only few results about systems initialized with random initial coloring. Given graph G We 
consider random initial set of seeds S p = S P (G) as the set containing any vertex of G with prob- 
ability p (independently on other vertices). Naturally, for every fixed graph G this gives us some 
probability that the random set of seeds is a dynamo. Determining these probabilities analytically 
is quite hard (if not impossible) and moreover, it can be done numerically with sufficient accuracy. 
Therefore, we shall try to obtain asymptotic results of the form: assuming 4-regular graph with n 
vertices, random set of seeds S* 0,12 is dynamo with high probability (w.h.p.). On the other hand, 
go. 10 j g no j. a dy nam0 (w.h.p.). 

In many situations, to determine the minimal value of p that S p is (w.h.p) dynamo is trivial. 
For example, assuming Erdos-Renyi random graph G(n,p') for fixed p', S p forms a dynamo almost 
surely (a.s.) if p > 1/2 + e (where e > is arbitrarily small constant) and does not form a dynamo 
a.s. if p < 1/2 — e ( this follows from Chernoff bounds and Markov inequality). If in the same 
model we allow p' to be dependent on n, but not to be significantly decreasing, the same results 
can be easily derived. Finally, if p' decreases significantly with n (for example p = c/n), the graph 
contains (with probability tending to some 7 > 0) some isolated vertices. Therefore until p is not 
quite close to 1 (so close that all isolated vertices are a.s. seeds) we can not say that S p is w.h.p. 
dynamo. Another trivial example is the toroidal mesh under strict majority contamination rule 
(that is the white vertex turns black only if it has at least three black neighbors). In this case, 
any square of size 2x2 forms an immune subgraph (that is such subgraph that turns black only 
if some of its vertices are seeds). Therefore, the sampling probability that would form dynamic 
monopoly must prevent all such squares from being colored entirely white. Once again, this can 
be done only for p very close to 1. 

However, the motivation presented in previous text leads to considering p to be a fault proba- 
bility of a node in the network (or something alike). Therefore we can assume 1 — p not to be very 
small. This makes the cases where p — > 1 non-realistic. 

In this work we shall examine dynamic monopolies with random initial condition on two types 
of underlying graphs. In the section 2 we focus on the toroidal mesh. The results are quite 



surprising - S p containing only a fraction o(l) of all the vertices a.s. form a dynamo. To be more 
concrete, we shall show that there exist constants a, (3 such that if p > a/ln(n), then S p a.s. is a 
dynamo. Similarly, when p < (3/ ln(n), then S p a.s. is not a dynamo. At the end of the section we 
present our attempt to solve the problem numerically and we show "measured" a and /3. As it was 
said earlier, to form a dynamo on random graph G(n,p') we need p to be close to 1/2. This makes 
the results about the toroidal mesh even more interesting. The natural question arises, whether 
the reason that such low values of p are needed to form a dynamo comes from specific topology of 
toroidal mesh, or whether it is just implied by the fact, that the degree of the vertices is constant 
and low. This is the motivation for the section 3, where we investigate the same questions but on 
random 4-regular graphs, that is the graphs with all the vertices having degree 4. We show that 
if p > 0.12, then the random initial coloring a.s. forms a dynamo. Similarly, if p < 0.10, then the 
random initial coloring a.s. does not form a dynamo. 

2 Toroidal mesh 

2.1 Preliminaries 

By toroidal mesh of size n we shall consider an undirected graph G = G(V, E) consisting of n 2 
vertices labeled as V[i,j] for < i, j < n. The set of edges consists of all pairs (V[i, j], V[(i + 1) 
mod n, j]) and (Vfi,^], V[i, (j + 1) mod n}) for all < i, j < n. By rectangle of size w x h located 
at position x,y we denote the subgraph of G induced by vertices V^j] for x < i < x + w and 
V < j < y + h. Such rectangle will be denoted as V[x,y,w,h]. The set of all vertices of the 
rectangle R having degree (in R) lower than 4 is called circumference of R. Rectangles with size 
w x 1 or 1 x h are called lines, rectangles with w = n or h = n are called vertical or horizontal 
stripes. Note that circumference of a stripe consists of two closed lines. 

Let ip : V — > N and S C V be the set of seeds. By coloring process induced by S we shall 
denote progression of sets Af(S, G, ip) = S , Si, . . . where S — S and Si C V represents the set 
of black vertices in the i-th step of the coloring process. The relation between Si and S^+i is that 
Si + \ = Si U B i+ i, where B i+ i C V is the set of all vertices v that are white in the i-th step of the 
coloring process (that is, v ^ Si) and they have at least ip(v) neighbors that are black in the i-th 
step of the coloring process (that is, they belong to Si). Furthermore, we shall set B = S and 
we shall say that vertices from Bi are turned black in time i. Vertex v is said to be turned black 
(by the coloring process) if it is turned black in some time. Since the object of our interest is the 
4-regular graph under the simple majority scenario, we shall put throughout the paper ip(v) = 2 for 
all v G V. If for some i, the coloring Si contains all vertices V, then we say that the corresponding 
S is dynamic monopoly or dynamo (with respect to G and ip). 

Vertex that belong (does not belong) to set of seeds S shall be called S- vertex (non-S*- vertex). 
Similarly, rectangle (line) consisting of S vertices (non-S'-vertices) is called .S-rectangle (S'-line, 
non-S 1 - rectangle, non- S'-line). 

We want to analyze the coloring process with random initial condition. For this purpose, we 
define the random set of seeds S p = S P (G). This is obtained as follows: Every vertex v of G is 
contained in S p with probability p, independently on the other vertices. All possible set of seeds 
form together with adequate probabilistic measure a probabilistic space of random set of seeds 
S p = S(G) P . We put S p G S p throughout the text. Naturally, probability p can depend on n. 
Since our aim is to obtain an asymptotic description of the behavior of the system, all limits we 
use are for the case n —¥ oo. We shall say that event A happens with high probability (w.h.p.), if 
linin^oo F(A) — > 1. Alternatively we shall say that A happens almost surely (a.s.). 

For abbreviation we shall put I = l(n) = |_ln(n)J and q = 1 — p. 



2.2 Lower bound 

In this section we show that unless p is high enough, there w.h.p. exists covering of G by cages 
(for the definition of cage see below), that prevents the black vertices from spreading. The result 
of this section is stated in the following lemma and proved at the end of the section. 

Theorem 2.1. If p < 0.012/ln(n) ; then the random set of seeds S p w.h.p. is not a dynamo on 
graph G. 

In this section we shall assume n to be divisible by 21. We shall show that this is only a technical 
complication at the end of this section. 

Let G be a toroidal mesh of size n x n and let S p G S p be random set of seeds. From this 
graph we can derive corresponding graph G' and corresponding set of seeds S' as follows: G' is a 
toroidal mesh of a size n/2 x n/2. We shall say that V'[x,y] corresponds to foursome of vertices 
V[2x, 2y],V[2x + 1, 2y],V[2x, 2y + 1], V[2x + 1, 2y + 1]. Vertex v' G V belongs to S' if at least 
one of its corresponding vertices is S^-vertex. Note that probability p' that the vertex in G' is 
S'- vertex is 

p' = 1 - (1 - p) A < Ap 

and for p — ¥ we get p' ~ Ap. 

Definition 2.2. By cage we denote the circumference of some rectangle R in G' consisting of 
non-S"-vertices only. In more detail, for any integers x, y, w < n, h < n if the vertex of G' on the 
coordinates [x + i,y],[x + i,y + h],[x,y + j],[x + w,y + j] is not member of S' for all < i < w and 
< j < h, then these vertices form cage in G' and vertices corresponding to these form cage in G. 
By the interior of the cage we denote those vertices of R that do not lie on the circumference. 

Lemma 2.3. If there exists a set /C of cages such that 

• every v' G S' p belongs to the interior of some cage K G /C; 

• for Ki, Kj G /C and Ki ^ Kj the interiors of the cages Ki and Kj are disjoint, 

then S p is not dynamic monopoly on G. 

Proof. We can assume all seeds to be "trapped" inside cages which they can not escape from. 
More rigorous proof can be done by induction on the number of steps of the coloring process. □ 

Let us consider rectangle R = V'[x,y,w, h] of G'. We shall call this rectangle good, if 

• w > h and there exists a non-S"-line V'[x, y + i, w, 1] for some < % < h of R, or 

• w < h and there exists a non-S"-line V'[x + i, y, 1, h] for some < % < w of R. 

If none of the above holds, we shall call R bad. 

Lemma 2.4. Let a be arbitrarily low positive number, (5 the root of the equation 

ln(l - e~ 12/3 ) = -2 

and let p be equal to ((3 — a)/ln(n). Then G' a.s. does not contain bad rectangle of size I x 31 or 
31 x /. 



Proof. As was mentioned before, p' ~ 4p = 4(/3 — a)/ln(n) The probability that the rectangle 
with width w and height h, w > h is bad equals to: 

Poncr C c = (l-(l-p'T) h . 

There are 2n 2 different rectangles of given dimensions, therefore the probability that there will be 
at least one bad rectangle of given dimensions can be upper-bounded using union bound by: 



^ 2 - -, 

J onc rcc 



Pall recs < 2n ■ p 

By letting w — 31, h — I, and using (1 — p') w ~ e- 12 (P~ a ) we immediately get 

Pall rccs ->■ as U ->■ (X). 

n 

Now we show, how it is possible to satisfy preconditions of Lemma 2.3. Assuming that n is 
divisible by 21 we divide G' to n/(2l) vertical stripes of width /: 

Z = \/'[0,0,/,n/2] 
Z 1 = V'[l-l,0,l,n/2] 
Z 2 = V'[2l-l,0,l,n/2] 



In each of these stripes we find closed path consisting of non S' vertices that goes "around" this 
strip. This will be done as follows: In every stripe Z^ we construct a sequence of rectangles 
Ri,ii Ri,2i ■ ■ ■ as in the top-left part of Figure 1. Due to Lemma 2.4 in each R it j there is a non-S"- 
line Li j. Suitable parts of these lines form desired path as is shown in the figure. Let us denote 
the path that goes around the stripe Zj by Pi and let Qi = {Lij\j G N}. Let us now consider all 
horizontal lines from Qi U Qi+\. These clearly divide the whole area between paths P\ and P 2 to 
cages with disjunct interiors. By repetition of this argument we can show that the area between 
any two paths Pi and P(i+\) mo d n can be covered by cages with respect to preconditions of Lemma 
2.3. For proof of Theorem 2.1 it is therefore sufficient to calculate appropriate /3 and p that would 
satisfy precondition of Lemma 2.4. 

The very last thing needed to finish the proof of Theorem 2.1 is to deal with such toroidal 
meshes that do not have its size divisible by 21. We shall use the following intuitive lemma: 

Lemma 2.5. Let H be the toroidal mesh, Sh the set of seeds, and let t be such a time after which 
no vertex changes its color to black. Let S* H be another set of seeds satisfying Sh C S h and let 
us denote D := S H — Sh- Then every vertex that turns black in J\f(S H ,H,2) after time u > t is 
distant at most 2u from some vertex of D. 

We assume that n = 2lk + r for some k and r < 21. As before, we divide G to vertical stripes 
with width 2/ (note that since we do not assume n even as before, we are working with G instead 
of G') and one special stripe with width r. As above, whole graph G up to the special stripe can 
be covered by cages. Every such a cage has the size of its interior limited to (2/) • (6/) = 12Z 2 . 
Therefore, if the special stripe contains only non-S^- vertices, then in time higher than 12Z 2 no 
vertex is turned black. Therefore the preconditions of the Lemma 2.5 are satisfied with H = G, 
D equal to the set of all ^-vertices from the special stripe, S H = S p , Sh = S H — D and t = 12/ 2 . 
Therefore we know that every vertex that turns black after time t must lie within the distance 
2t from some black vertex from the special stripe. Note that dividing G to the stripes can be 



Figure 1: Left: Three consecutive stripes. In every (dashed) rectangle there exists a non S' line 
(bold lines). These lines form path around the stripe. Right-top: Black rectangle R will expand 
itself to ABCD because of incrementing vertices V{. Right-middle: Every square 2/ x 21 (dashed) 
w.h.p. contains this configuration of non S' lines (bold lines). Bottom: results of a numerical 
study. 
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done in many ways - the special stripe can be located at arbitrary position. If we assume that 
special stripe equals to l / [0,0,r, n], then every vertex that turns black after time t must has its 
x coordinate within the interval (—2t,2t + r). On the other hand, if the special stripe equals to 
V[|ji/2J, 0, r, n], then every vertex that turns black after time t must has its x coordinate within 
the interval {\n/2\ — 2t, \ji/2\ +2t + r). Since t = o(n) and r = o(n), these two intervals do not 
overlap and therefore no vertex turns black in time t + 1. This implies that no vertex turns black 
anymore and therefore the initial coloring S p is not a dynamo. 

2.3 Upper bound 

The result of this section is stated in the following theorem: 

Theorem 2.6. If p > 1.65/ ln(n), then the random set of seeds S p w.h.p. is a dynamo on graph 
G. 

We want to illustrate the main idea of the proof first. The rigorous proof is presented after few 
technical lemmas. Let R be such square that in some step of the coloring process consists of black 
vertices only. If in each one of the four dashed lines (Figure 1 top-right) there exists at least one 
seed (as are Vi,v 2,^3,^4 in the figure), then also the square ABCD will be black at some time. 
We shall call vertices t> 1 ,t>2,t>3,t>4 incrementing with respect to R. Expansion of R caused by 
incrementing vertices can possibly continue until whole G turns black. We shall show that we can 
choose R such that with probability tending to one this really will happen. 

Let h = h(n) such that h(n) is odd. Let us divide G to the set of \n/h\ 2 disjoint squares Ui 
such that the size of Ui is h x h. The event that G turns black at some point in the coloring process 
is implied by existence of specific square Uk that satisfies: 

• Ti: A square consisting of single vertex v in the middle of Uk can grow due to the existence 
of incrementing vertices from S p in such a way that in some time of the coloring process the 
whole Uk will be black. Note that the square of size lxl can grow due to incrementing 
vertices to the square of size 3x3 even if it is not black. 

• T2'. Square Uk will grow due to the existence of incrementing vertices from S p , and at some 
time whole G turns black. 

Let us note that although it is tempting to combine the conditions r\ and T2, it can not be 
done easily. The reason is that for two different U,Uj, the corresponding IVs are independent 
events. On the other hand, the "spreading" of the square C/j to the whole G is influenced by (not) 
spreading of other squares Uj. 

We shall estimate probabilities of conditions Ti and T2 right after stating following technical 
lemma. 

Lemma 2.7. Let p > c/l for some constant c. Then every line of length L := I 3 contains w.h.p. 
at least one S p -vertex. 

Proof. The probability that the fixed line of length L is non-S^-line is 

Pone line = <? ~ e F = e 

With the use of the union bound we get that the probability that at least one line of length L is 
non-SMine can be bounded as 

Pall lines _ ^^ ' Pone line ' "■ 

□ 



Lemma 2.8. Let us assume that h — > oo. Then the condition 

q ■ i\ 2 

ln(n) ln(l/p) — ln(/i) — >■ oo 

op 

w.h.p. implies that I\ happens for arbitrarily large number of Uk 's. 

Proof. Let Aj denote the event that for given C/j the condition Tx holds. Let X be the random 
variable that equals to the number of Aj's holding. There are at least \n/h\ 2 ~ (n/h) 2 different 
C/j's and therefore the expected value of X satisfies the relation E(X) ~ (n/h) 2 F(Ai). Since the 
AiS are independent, E(X) — )■ oo implies that w.h.p. at least one of the A^s holds. Let us estimate 
¥(Ai) (for fixed U { ) as follows: 

P(A)>Kdd where Kdd = J] (1 - q 1 ) 4 (1) 

j>l,i is odd 

Furthermore we define 

Veven = Q (1 ~ <f ) 4 and Kll = V even ■ Vodd 

i>l,i is even 

Note that V odd < V cvcn but V^ dd > V^n ■ (1 - g) 4 = V cvcn ■ p 4 . Therefore V odd > V / Ku • P 2 - 
Further, we compute: 

oo oo oo 

Kn = n( 1 -^) 4 = ex p( 4 E ln ( 1 -^)) =ex p(- 4 E^+^ 2 v 2 +? 3 V3 

i=l i=l i=l 

oo oo oo oo 

= exp (-4 J2 E Q ki / k ) = ex P (" 4 E E «*/*' 

i=l fc=l fc=l i=l 

00 i A; °° i i 

V ^k' l-a k ) ~ exp v ~ 2-^ jfc " a -fe_ i, 



exp | 



fc=i ^ fe=i 



where in the second step we used Taylor's expansion for ln(l + x). Furthermore, we use Bernoulli's 
inequality: 

q - k -\ = (\ + z ) k -\> kz 
where we put 1 + z — 1/q. Using this we obtain 



oo oo 

mw/4 > - E i = -i E p = -' 2 /6 ■ i = -* 2 /e 

fc=i fc=i 

Finally E(X) — >■ oo can be written as: 



(n/h) 2 -F(A t ) > (n/h) 2 Jexp(-A- 7 ^-') ■ p 2 ^ oo 

Calculating logarithm of the last relation proves the lemma. □ 

Proof. (Theorem 2.6). Let us choose h = 2/ 3 + 1 and let p be as in the statement of Theorem 
2.6. The direct calculation shows that the preconditions of Lemma 2.8 are satisfied, and therefore 
there a.s. exists Uk that satisfies the condition Ti. Moreover, since h > I 3 , Lemma 2.7 implies that 
there w.h.p. exist all incrementing vertices needed for the condition T 2 to hold. □ 



Let us note that both Theorems 2.6 and 2.1 can be sharpened a little. When improving Theorem 
2.1 we can show (with some additional effort) that for p < 0.018/ ln(n) every square of size 11 x 21 
a.s. contains four non-S^-lines such as in the Figure 1 (middle-right). Existence of such structures 
is sufficient for showing the existence of covering of G by cages as required by Lemma 2.3. For 
improvement of Theorem 2.6, let us note that we do not actually need all incrementing vertices 
whose existence is included in (1) in Lemma 2.8. For successful "spreading" of black vertices it is 
sufficient if only approximately half of them exist. Let us recall that in Lemma 2.8 we started with 
lxl square and this square potentially grew until the whole Sk became black. The expansion of 
this square of size % x i to (i + 2) x (i + 2) requires the existence of the set of four incrementing 
vertices Ui. Another set of incrementing vertices Ui then causes expansion to size (i + 4) x (i + 4). 
But in almost all situations, XJi can cause expansion from % x % to (i + 4) x (i + 4) directly - without 
need of U±. This leads to p > 0.83/ ln(p). 

This considerations sharpen the results, but it is still an interesting question, whether there 
exists a threshold t such that if t' > t and p > t'/\n(n), then the initial random set of seeds is 
a dynamo and vice-versa. Sadly, we are not able to prove the (non)existence of such threshold 
and determine its possible value. Therefore we present a numerical study that tries to answer 
these questions at least partially. By p z (n) let us denote the probability p, such that the random 
set of seeds S p is a dynamo with probability z. Results of the numerical study are presented in 
the Figure 1 (bottom). The circles present values of po.s(^)- The top and the bottom points of 
the "error bars" correspond to values Po.ob(n) and Po.m{n). Size of 90% confidence intervals for 
measured values is comparable with the diameter of the circles. 

We conclude that threshold t is likely to exist and its value is near 0.25. Proving this statement 
rigorously remains an open problem. 

3 Dynamos on random 4-regular graph 

In the previous section we analyzed a particular 4-regular graph. In this section we try to solve a 
similar task for a random 4-regular graph under simple majority scenario. 

Random regular graphs are very well studied. Pioneering studies in this field were brought by 
Bender and Canfield [12], Bollobas [13] and Wormald [14, 15]. A systematic research in this area 
grew enormously since then, partly driven by applications in many areas such as computer science. 
For exhaustive survey see [16]. To make this paper self-contained, we now present the definition 
and basic properties that are necessary for further reading. 

We use Q n< d to denote the uniform probability space of ci-regular graphs on the n vertices 
{1,2,..., n} (where dn is even). Sampling from Q n & is therefore equivalent to taking such a graph 
uniformly at random (u.a.r.). Another possible way how to define (the same) probabilistic space is 
this: We construct graph with n vertices (named as before) and no edges. Then, for every vertex 
i we construct d slots v^i, . . . ,Vi t d- A perfect matching of these slots into dn/2 pairs is called a 
pairing. A pairing P corresponds to a multigraph (with loops and multiple edges permitted), in 
which two vertices i and j are connected, if there exist two different slots Vi yX and v^ y that form a 
pair from P. Although this process can lead to graph that is not simple, it is quite easy to show 
that the probabilities of obtaining two simple graphs G\ and G<i are equal. Therefore if we reject 
every graph that is not simple and repeat the whole process until simple graph is found, we obtain 
the same probabilistic space Q n ^ as before. 

A pairing can be selected u.a.r. in many different ways. In particular, the slots in the pairs 
can be chosen sequentially. At any stage, the first slot in the next random chosen pair can be 
selected using any rule whatsoever, as long as the second slot in that pair is chosen u.a.r. from 
the remaining slots. For example, one can insist that the next chosen slot is the next one available 



in any pre-specified ordering of the slots, or comes from a cell containing one of the slots in the 
previous chosen pair (if any such slots are still unpaired). We shall call this the independence 
property of the pairing model. 

In what follows G G Q n ,A will be random 4-regular graph and I = \]n(n)\ as before. Also, we 
consider the same coloring process as the one described in preliminaries to section 2. 

Let us note that we do not need to restrict our interest to simple graphs. The coloring process 
is well-defined also for multi-graphs. Bender and Canfield showed that the probability of obtaining 
a simple graph by the pairing-slot process is asymptotically exp((l — d 2 )/4) {d being the degree of 
a vertex), which is a value close to 2% for 4-regular graph and high number of vertices. Therefore 
if we prove that S p is not a dynamo w.h.p. on random 4-regular multi-graph, it implies that S p is 
not w.h.p. dynamo on random 4-regular graph. 

In the future text, let G be random 4-regular (simple) graph of size n and H be random 
4-regular (not necessarily simple) graph of size n. 

Theorem 3.1. Let p > 0.12. Then S p w.h.p. is a dynamo on G. 

Proof. As was justified above, we shall work with H instead of G. We use some properties of 
the random regular graph that are summarized in [16]. The graph has only small number of 
short cycles. Therefore, the neighborhood of any vertex looks like a part of an infinite tree. For 
v G V(H), let T v be a rooted tree that is obtained as follows: Vertex set of T is a subset of V(H), 
vertex v is root of T. Tree T consists of h = |_2 • lg log(n)J levels. Vertex v forms 0-th level of the 
tree. The (i + l)-th level consists of all vertices that are neighbors (in H) to some vertex in the i-th 
level and that are not already contained in the j < i-th level for any j. If u, v are two vertices of 
i-th and (i + l)-th level respectively, then there exists an edge between u and v in T if there exists 
an edge between u and v in H and furthermore, there is no other vertex w than u from i-th level 
of T that is adjacent to v (in H). There are no edges within one level of T v . Note that although 
some edges from H are not present in T v , with the use of the independence property we can easily 
calculate that the probability that some vertex in i-th level (i < h) has degree lower than 4 is o(l). 
It is clear that if we run the coloring process on T v instead of H (but we still set the threshold 
ip to 2 for all vertices), then the probability that the root of T v turns black in coloring process 
induced by S P (T V ) is lower than the probability that the same vertex v turns black in H (in the 
coloring process induced by S P (H)). By Y(i) let us denote the probability that the vertex of the 
i-th level of T v turns black (in the coloring process induced by S P (T V )). We shall start analyzing 
the bottom level of the tree and proceed to the top. It obviously holds that Y(h) = p. For the 
(/ — i)-th level (i > 0) the situation is more complicated. In order to obtain some estimate for 
Y(l — i) we shall assume that the color of some vertex u from (/ — i)-th level can be changed only 
because of (at least) two black neighbors of u belonging to (I — i + l)-th level. This gives us the 
following estimation: 

Y(i-l)>p+(l-p).f 3 ,> 2 (Y(i)) (2) 



where 



f n ,=k(p) = ( ?) • P k ■ (1 - PT h and / n> > fc (p) = J2 /> 
^ ' i=k 



Numerical calculation shows that given p > 0.12 the condition Y(h — 100) > 0.999 holds. Here 
we change our approach and analyze higher levels analytically. We easily prove that for any p, 
q = 1 — p, it holds that: 

l-/ 3 »<5g 3 (3) 
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Let us denote £j = 1 — Y(i). Then we can rewrite (2) using (3) to the form: 

Y(i-l) >p+(l-p)-(l-5-e 3 t ) 

And finally, we get: e^\ < 5ef. Last inequality together with numerically obtained result e^-ioo < 
0.001 implies, that 

log(^-i) < log(5) + 3 log(ei) < 2 logfo). 

Therefore 

log(e ) < -(2 h ~ e W) = -log(n) 2 /e(l) < -21og(n) 

for sufficiently big n and finally Eq < n~ 2 . This means that the expected value of number of vertices 
from H that do not turn black is o(l). □ 

Theorem 3.2. Let p < 0.10. Then S p w.h.p. is not a dynamo on G. 

Before we present the actual proof of this theorem, let us recall several facts concerning the Galton- 
Watson branching process. 

Let Z be a probability distribution on the nonnegative integers. The Galton-Watson branching 
process with offspring distribution Z is (loosely) defined as follows. The 0-th generation consists of 
N particles. For t > 0, the (t + l)-th generation consists of the children of all particles in the t-th 
generation. For each particle, the number of its children has distribution Z, and is independent of 
all the other particles. 

Let X be a (Galton-Watson) branching process with N particles in the beginning, distribution 
Z and we shall denote the number of particles in generation t by X t and the average number of 
children of a particle by A (= E(Z)). Moreover, let 



oc 



N 



1-A 



i=0 

Lemma 3.3. For all t > 0, 

E(X t ) = N ■ A* and if A < 1 then E(X sum 

Furthermore, if X < 1, and a 2 is the finite variance of Z , then 

P(X sum >2-E(X sum )) = o(l). 

Proof. The proof of the first and the second statement is trivial. For the proof of the last statement, 
we use the first part of the lemma and variance estimation var(X sum ) < Na 2 /(1 — A) 2 (see for 
instance [11]). Therefore Chebyshev inequality is applicable, leading to the required result. □ 

Now we are ready to prove the main result of this section. 

Proof. (Theorem 3.2) As was justified above, we shall work with H instead of G. Although an 
approach quite different from that used in previous section is to be used for proving the Theorem 
3.1, there is one similarity: we first numerically compute how the coloring process evolves (w.h.p.) 
up to the certain time and only then we solve the problem analytically by finding analogy with 
easily solvable branching process. 

In Table 3 we present a randomized algorithm that simulates the coloring process on random 
4-regular multi-graph. The algorithm returns 1 with the same probability as is the probability that 
S p is a dynamo on random 4-regular multi-graph. Let us recall that the analysis of the coloring 
process on multi-graphs is sufficient for proving the statement of this theorem. 
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Table 1: Sampling-Coloring (SC) Algorithm 



1. 


Sample the set B and let _R := V{H) — B . 




2. 


Sample the neighbors of all vertices from B . Let Mi be those vertices from R 




that have one neighbor in Bo and let B\ be those vertices that have two or 


more 




neighbors. Let R\ := Ro — B\ — M\. 




3. 


Let % :— 1. 




4. 


Sample the neighbors of B,i and denote by T the set of them. 




5. 


Let Tl C Tj fl i?j be the vertices that have at least two neighbors in Bi. 




6. 


Let Tj 2 := T D M*. 




7. 


LetT 4 3 :=T-T/-T 2 . 




8. 


Let B i+1 := T} U T 2 . 




9. 


LetM m :=(M l -T?)UT i 3 . 




10. 


Let i? i+ i := i?j — £> i+ i — M i+1 . 




11. 


Let % :— % + 1. 




12. 


If Bj = and i^ = return 1; if T^ = and i^ ^ return 0. 




13. 


Go to line 4. 





The main idea behind the algorithm is that the edges in H are not sampled at once. Instead, 
in the i-th step of the coloring process we sample only edges corresponding to slots that belong 
to black vertices and have not yet been sampled (that is, Bi). By doing this, the sampling of 
the graph H takes place "in the same time" as the black color spreads. Because of independence 
property, this leads to the same probabilistic space of coloring processes as if we first sample the 
whole graph H, then the random initial coloring and only then start with the coloring process. 

The meaning of used symbols is: Bi (as defined in the preliminaries) is the set of vertices that 
turns black in the i-th step of the coloring process; Mi is the set of white vertices that have exactly 
one black neighbor in the i-th step of the coloring process (excluding Bi); Tj are neighbors of B{ 
and they are split up to the sets Tf (k = 1, 2, 3) according to their future contribution to one of 
following sets: B i+1 , B i+1 , M i+1 . 

We now calculate some estimations for the cardinalities of Bi, Ri and Tf. For abbreviation we 
shall use x for \X\ for all used sets. We shall restrict the considerations to % < 60 and p = 0.10. 
Note, that under these assumptions all cardinalities are O(n). We define Zi as the ratio of number 
of empty slots of the vertices from Bi to the total number of empty slots in the i-th step, that is: 



Zi < 



2-bi 



4-r. 
for i > and z = \B \/\V{H)\. By f n> > k , f n 
of Theorem 2.6. Now the following estimations hold: 



(4) 



h 3 • mi + 2 • bi 

& we shall mean the functions defined in the proof 



Lemma 3.4. For abbreviation, let X := t 
E(X)=r fc / 4) >2(* fc )(l + o(l)) 



\. Then, if r k ,bk = 0(n), then 
and var(X) < E(X)Q(l) + E 2 (X)o(l). 



Proof. Let X = ^ JQ where X t is an indicator random variable indicating that i-th vertex from 
Rk has at least two neighbors from B k , that is, it belongs to T k x (and consequently to -Bfc+i)- The 
probability that this happens (that is, Xi = 1) is clearly /4,>2(zfc)(l + o(l)). There are r k such 
vertices, which proves the first part of the lemma. For the second part, let us compute: 



var(X) = E(X 2 ) - E 2 (X) = J^E(XjX j ) - E 2 (X) = J^E(X i X i ) + E(X) - E 2 (X). 



1,3 



i¥=j 



12 



Clearly for i ± j holds E(X i X j ) = (fa,> 2 (z k )) 2 (l + o(l)). There are r k (r k - 1) 
sum. The second part of the lemma follows by straightforward calculation. 



such terms in the 

□ 



From this lemma and from Chebyshev inequality we have, that if r k ,b k = Q(n), then with 
high probability t\ = r k f^> 2 {. z k)l, where 7 is some constant that is arbitrarily close to 1. Quite 
similarly, we can derive such relations for cardinalities of all sets of SC algorithm. This gives us: 



h -- 


= P- 


n ■ 7 


r = 


= n - 


-b 


&i = 


= r 


■ fa,>2(z ) ■ 7 


ni = 


= r 


■ fa,=i(z ) -7 


7"i = 


= r 


— bi — mi 


t] - 


= n ■ 


■ U,>2( z i) ■ 7 



'■1 

t 3 

'■1 

b i+ i 

rrii+i 

r i+ i 



mf fs,>i(zi) -7 
mi- f 3 ,>i(zi) -7 

Ti- fi,=l{Zi) -7 

*} + *? 

rrii + t 3 - t[ 2 

ri - b i+ i - (m i+ i - mi) 



For proving that the coloring does not form a dynamo we need to upper bound rrii and bi and lower 
bound rj. Therefore we put Zi as big as possible (see (4)), and furthermore we put 7 = 1.0001 
in all equations above except for the case of tf, where we use 7 = 0.9999. With the use of this 
approach we obtain the following results: 



60 

m 60 < 0.36n b m < 10" 8 n r 60 > 0.41n I f] B J < 0.24n 

i=0 

The statement of the theorem looks now very persuasive - the fraction of vertices that turn black 
in the 60-th step of the coloring process is lower than 10~ 8 . However, we can not continue with 
evaluation of the above equations indefinitely, since we need all quantities to be at least @(n). 
Otherwise, it would be hard to bound 7 present in these equations. However, we can use the 
analogy with the branching process discussed by Lemma 3.3. Let us now consider the following 
conditions: 



m k > 0.37n b k > 10~ 6 n r k < OAOn I p| b\ > 0.25n 



(5) 



We prove that w.h.p. there is no such k that would satisfy a condition from (5). We already 
concluded that this is true for k < 60. Let us fix some coloring process on the graph H such that 
at least one condition from (5) holds and let k' be the minimal k for which this happens. It is 
easy to calculate that until no condition from (5) holds, the expected number of vertices that turn 
black in the (60 + fc)-th step (k > 0) can be upper bounded by the expected number of particles 
that exists in the k-th generation of branching process with parameters N = 10~ 8 n and A = 0.9. 
This gives us the equality E(X, um ) = 10~ 7 . Note also that any of the conditions from (5) implies 
that the number of the particles in the branching process exceeds 100 iV = 10 • E(X sum ), which by 
Lemma 3.3 happens with probability o(l). Therefore, w.h.p. less than 40% from all vertices turn 
black. □ 

Finally, there are rounding errors in our numerical analysis. Influence of these errors can be 
estimated by multiplication (or division) of the result of each numerical operation by some number 
close to 1. In the proof of Theorem 3.2 this is clearly hidden in our choice of 7; in the proof 
of Theorem 3.1 we can simple multiply the right side of (2) by 1 — 10~ 6 which produces more 
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significant "error" than the error of the machine. However, the result Y(h — 100) > 0.999 still 
holds. 
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